Perturbed breakup of gas bubbles in water: Memory, gas flow, and coalescence 
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The pinch-off of an air bubble from an underwater nozzle ends in a singularity with a remark- 
able sensitivity to a variety of perturbations. I report on experiments that break both the axial 
(i.e., vertical) and azimuthal symmetry of the singularity formation. The density of the inner gas 
influences the axial asymmetry of the neck near pinch-off. For denser gases, flow through the neck 
late in collapse changes the pinch-off dynamics. Gas density is also implicated in the formation of 
satellite bubbles. The azimuthal shape oscillations described by Schmidt et at, can be initiated 
by anisotropic boundary conditions in the liquid as well as with an asymmetric nozzle shape. I 
measure the n — 3 oscillatory mode, and observe the nonlinear, highly three-dimensional outcomes 
of pinch-off with large azimuthal perturbations. These are consistent with prior theory. 
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I. INTRODUCTION 

The breakup of an underwater bubble, whether in a 
swimming pool, crashing surf, or even in the processes of 
a paper mill, exemplifies the rich underpinnings of even 
the most fleeting phenomena. In a brief moment, a thick 
neck of air collapses to a microscopic radius and breaks 
— a process that can be modeled as the formation of a 
singularity, where physical quantities, including pressure 
and velocity, diverge [H-ll]- It was previously thought 




FIG. I: Air bubble pinch-off. The bubble is blown from a 
6 mm-diameter circular hole cut in a plate that was carefully 
leveled. The plate reflects the image of the bubble. The bright 
spot in the dark bubble is due to back-lighting. 
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that all breakup singularities, from dripping faucets to 
underwater bubbles, are universal, so that as the pinch- 
ing neck of fluid approaches the singularity, its shape 
and dynamics are determined solely by the two fluids' 
material parameters, and initial conditions are forgot- 
ten. Recently, it was discovered that the singularities at 
the end of certain pinch-offs instead retain information 
about their beginnings [!, 0] . 

In the case of an air bubble pinching off from an un- 
derwater nozzle, illustrated by Fig. [TJ the singularity has 
a detailed memory of the bubble shape as pinch-off is 
beginning, which may change the character of the singu- 
larity formation, or even disrupt it entirely @, Q ■ These 
developments raise two central questions: What informa- 
tion is preserved at the singularity? What terms in the 
dynamics can perturbatively break the symmetry of the 
singularity and thus disrupt the singularity formation? 
Studies of such a rich phenomenon as air bubble pinch-off 
can thus provide new insights into dynamical singulari- 
ties, and the free-surface flows that develop them. 

This paper reports on experiments that investigate fur- 
ther what information is preserved at the singularity. It 
focuses on the way the azimuthal and axial symmetry 
of pinch-off are broken by initial and boundary condi- 
tions [1,0], and by the gas inside the bubble [9kTTT | . 

Past and present studies of perturbed breakup begin 
with models of the fully symmetric pinch-off of a void in 
an inertial fluid, and the experiments that address these 
models. Previous investigations of this kind of pinch-off 
have focused on the shape of the cavity near the sin- 
gularity, and the extent to which the evolution of the 
neck differs from that of a universal self-similar singu- 
larity. These studies concentrate on determining the ra- 
dius -R m in of the narrowest part of the neck, measured 
as a function of the time r = t* — t remaining until 
the time of disconnection t*. The model of Longuet- 
Higgins et al. uses the approximation that conditions 
vary slowly with vertical position z, so that the pinch- 
off dynamics can be approximated in a two-dimensional 
horizontal slice m . It shows that the inertia of the water 
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dominates the asymptotic dynamics, so that surface ten- 
sion, viscosity, and gas may be neglected, and hydrostatic 
and Bernoulli pressures drive the collapse. In this limit, 
the evolution should approximately follow the power law 
-Rmin oc T 1 / 2 , a finding borne out by simulation and ex- 
periment nmm. 

Subsequent theoretical and simulational work has 
shown that when the three-dimensional shape of the ax- 
isymmetric neck is considered, the exponent of this power 
law may evolve during the collapse, varying with the log- 
arithm of T Such changes over many orders of 
magnitude cannot be measured in experiment, but these 
models are consistent with experiments that have mea- 
sured a single value for the exponent that is slightly 
higher than 0.50, typically 0.56 % The present 
work finds similar values for the exponent. 

In this paper, I first discuss the methods used to ob- 
serve and measure pinch-off, to vary the composition of 
the gas, and to break the azimuthal symmetry of the 
bubble. This is followed in Sec. IHII bv results on the role 
of gas in altering the neck shape and creating satellite 
bubbles. In Sec. IIV1 I report on experiments in which 
the azimuthal symmetry of the neck is broken by initial 
or boundary conditions. I observe vibrations of the neck 
shape, and in some cases, new shapes and topologies. 



II. EXPERIMENTAL METHODS 
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My experiment in its simplest form is represented in 
Fig. [2ji. Gas is blown slowly from an underwater nozzle, 
adding to the volume of a bubble until surface tension 
can no longer hold it at the nozzle, so that it begins to 
pinch off y . Gas flows from a pressure-regulated source 
through a needle valve, and so is added at a constant rate 
to the combined volume of the bubble and the tubing 
connected to the nozzle (~9 mL). The top of the nozzle 
is fitted with one of a series of interchangeable aluminum 
plates, each with an orifice of a specific size and shape 
cut out of it, from which the bubble emerges. 

Prior experiments have shown that underwater bubble 
pinch-off has a remarkable sensitivity to any tilt of the 
nozzle orifice from the horizontal @ . More generally, the 
shape of pinch-off is a function of the contact line at the 
base of the bubble, where the gas- water interface termi- 
nates. Therefore the position of this contact line must be 
controlled carefully. First, the line's position with respect 
to the nozzle must be fixed. To this end, the inside of the 
nozzle is machined to reduce its roughness, and made hy- 
drophobic with a coating of PTFE ("Fluoroglide," Saint- 
Gobain); the PTFE is sanded off the plate's top surface to 
make it hydrophilic. This design ensures that the bub- 
ble interface stays pinned to the top rim of the nozzle 
orifice. Second, the orientation of the nozzle itself must 
be controlled, so that the top of the nozzle is in a plane 
perpendicular to gravity. A leveling adjustment is there- 
fore performed at the beginning of each experiment, by 
changing the inclination of the water tank in increments 



FIG. 2: (a) Schematic top view of the experimental appara- 
tus. Bubbles are blown slowly from a nozzle, back-lit, and 
imaged synchronously with two high-speed cameras. Several 
nozzles of different shapes and sizes may be fitted to the ap- 
paratus, (b) Variant of the gas supply of (a) for generating 
bursts. Gas is slowly released through a needle valve until the 
bubble reaches a preset state; a pressure sensor then triggers 
the release of the gas held in the reservoir. 



as small as 0.01° until pinch-off has bilateral symmetry 
when viewed from both the front and side. Because of 
these measures, azimuthal asymmetry of pinch-off in my 
experiments is due exclusively to the manufactured shape 
of the nozzle orifice, or to deliberate tilting. 

Bubbles are back-lit and imaged with high-speed video. 
The use of two cameras, positioned at a right angle 
to each other, permits a measurement of the neck's 
azimuthal asymmetry. Most experiments reported in 
this work were imaged with high-speed video cameras 
("Phantom" vl2.1 or v7.3 from Vision Research, Inc.) 
at up to 183,000 frames per second with a 1 /is expo- 
sure. A long-distance microscope ("InfiniMax" by In- 
finity Photo-Optical) on each camera typically provided 
a magnification of 3.8 /^m/pixel. Movies were lit with 
low-infrared continuous tungsten lights ("DedoCool" by 
Dedotec), which prevent significant heating of the wa- 
ter. Some single exposures were lit with a 300-ns arc 
flash ("Palflash" from Cooke Corp.), triggered when the 
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disconnecting bubble interrupted a laser beam. 

Each resulting movie is analyzed by first selecting a 
threshold pixel value halfway between the most common 
black and white pixel values in the movie, then using 
this threshold to find the positions of the neck edges for 
each row of pixels, with linear interpolation between ad- 
jacent pixels. Pixels where dirt or a camera malfunc- 
tion would cause erroneous results are excluded from the 
analysis. This process yields the neck profile, which is 
the half- width R as a function of height z. Random er- 
ror in the edge position, due to both image noise and the 
edge-finding process, is manifest as fluctuations in the 
computed midline of the neck, which would ideally be a 
straight or gently curving vertical line. These fluctua- 
tions were found to be very small, typically well below 
0.1 pixels. To find i? m ; n from R(z), a parabola is fitted 
to a small region around the minimum. Distance scales 
were determined by imaging a reference cylinder, 2 mm 
in diameter, at the focal plane. 

To analyze scaling behavior of the minimum neck ra- 
dius -R m i n with time, the R m i n (t) data from each individ- 
ual movie are fitted with the function: 

Rmin=P(U-t) a . (1) 

To determine the range of exponents a that is consistent 
with a given movie, a number of trial fits are performed 
by varying t* as an independent parameter, and observ- 
ing the quality of each resulting fit and the value of a that 
it yields. When judging the quality of a fit, I assume a 
systematic error for i? m i n of ±0.5 pixels (as opposed to 
the much smaller random error described above). This 
reflects the possibility that the edge-tracing algorithm, 
combined with the optical system, might systematically 
misjudge the position of the edge. The resulting best 
value of a, and range of i? m in to include in the fitting 
process, are then used to determine t* and /3 for each 
movie in the data set, so that all data may be plotted on 
a single time axis, r = t* — t. 

A. Repeatability of quasistatic experiments 

My normal method of bubble inflation, used in the 
majority of experiments described here, involves slowly 
inflating the bubble until its buoyancy overcomes surface 
tension and it pinches off. When bubbles are inflated in 
such a quasistatic way, and the position of the nozzle and 
gas- water interface is well-controlled as described above, 
pinch-off is initiated with a precise bubble size and shape, 
so that these quasistatic experiments are remarkably re- 
peatable and reproducible. Nearly all of the plotted data 
sets in the present work are taken from at least two sepa- 
rate pinch-off events, and the experimental and analysis 
protocols described here were developed with the goal of 
ensuring consistency between results obtained from dif- 
ferent runs. 

Two further sources of non-repeatability that are ad- 
dressed by my experimental protocols bear mentioning. 



First, my most sensitive measurements of the asymmetry 
introduced by a tilted nozzle suggest that the water cur- 
rents excited by repeated quasistatic pinch-offs can have 
a minute, non-reproducible, systematic effect on this tilt- 
ing asymmetry at pinch-off. Waiting at least 90 s between 
bubbles eliminates this effect. However, these most sensi- 
tive experiments on tilting were the only type of measure- 
ment where I could find a systematic effect from faster 
quasistatic bubble production. For all other types of ex- 
periments presented in this paper, I was unable to iden- 
tify any systematic difference between 90 s/bubble and 
~10 s/bubble, although a small increase in random vari- 
ation between bubbles was observed at the faster rate. 
No data reported in this paper were collected faster than 
-10 s/bubble. 

Second, external vibrations affect the timing and sym- 
metry of pinch-off, and in some cases, seriously increase 
the random scatter of resulting data. Experiments were 
therefore conducted on a heavy optical table supporting 
as few moving components as possible, and with limited 
human or mechanical activity in the laboratory. The 
optical table's isolating suspension system could not be 
used, however, because it introduced error into the ta- 
ble's orientation with respect to the horizontal. 



B. Measuring asymmetry with 2 cameras 

Measuring azimuthal asymmetry with 2 cameras, as 
in Fig. [2j involves finding the difference between the 
cameras' simultaneous measurements of the neck profile 
R(z). This subtraction requires transforming the two- 
dimensional profiles, originally in units of pixels, to a sin- 
gle three-dimensional laboratory coordinate system with 
units of /im. Error in each camera's conversion from 
pixels to ^m (i.e., the magnifications) is the dominant 
source of systematic error in this transformation. (By 
contrast, the aforementioned ±0.5 pixel systematic er- 
ror in locating the true positions of the neck edges is 
common to the measurements of both cameras, and thus 
cancels out in the subtraction.) While the magnification 
of each camera may be measured by imaging a refer- 
ence object, this cannot be done with enough precision 
to then measure an asymmetry much smaller than the 
neck size. For example, in order to determine that a 
neck with i? m i n = 100 /im is 1 /im wider in one camera 
than in the other, each magnification should be known 
to better than 1 part in 100. 

Fortunately, it is possible to measure instead the ratio 
of the 2 magnifications with sufficient precision. The ra- 
tio may be worked out by using movie frames from each 
camera to measure the separation between the 2 cusps of 
air (connected to the bubble and to the nozzle) that re- 
main after breakup, as they retract vertically away from 
the pinch-off point. When repeated over many pinch-off 
events, this method is equivalent to measuring many ref- 
erence objects spanning a range of length scales. Further- 
more, unlike when using a manually positioned reference 
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object, the calibration and experimental measurements 
are necessarily in the same focal plane. Each camera's 
line of sight is kept sufficiently horizontal so that the 
mismatch between horizontal and vertical magnifications 
due to projection effects is negligible. This technique also 
yields the elevation of one camera relative to the other, 
which is required to compare data from the 2 cameras 
on a single z axis. Because this calibration process is ap- 
plied each time the optical system is adjusted, I am able 
to measure reproducibly asymmetries in the neck shape 
to better than 1 part in 500, as estimated by comparing 
similar measurements taken on multiple days. 



C. Control of gas 

Most of the experiments reported here used nitrogen as 
the gas. To control the gas density and the amount of dis- 
solved gas, some experiments were also performed with 
He and SFg gas in a glass vacuum cell that allowed the 
observation of pinch-off through front and side windows. 
A vacuum pump with a cold trap was used to lower the 
pressure to ~10 kPa, though observations were made only 
at ^20 kPa and above, because of the difficulty of main- 
taining a stationary contact line at the nozzle when gas 
density is low (see Sec. Qll] below). Water was degassed 
in a separate vessel before the experiment. It was further 
degassed in the experimental cell by lowering the pres- 
sure to ~10 kPa, then returning to atmospheric pressure 
by slowly bubbling the less-soluble He or SF 6 through the 
nozzle. This procedure was repeated at least once before 
a final pump-down to ^20 kPa. The nozzle was then 
leveled according to the procedure described above, after 
which recorded observations began. Throughout the ex- 
periment, pressure was raised solely by introducing gas 
at the nozzle. This ensured that the quantity of dissolved 
gas would increase minimally during the experiment, as 
SF 6 and He have low solubility in water (0.007 and 0.009 
by volume, respectively, as compared to 0.03 for air plj). 

D. Bursting apparatus 

To study the pinch-off singularity's memory of per- 
turbations, a nozzle with a non-circular orifice is used 
to break the azimuthal symmetry of the bubble shape. 
Most such experiments are carried out with quasistatic 
inflation to maximize reproducibility (see Refs. @, 0| and 
the present work). However, in the course of quasistatic 
inflation, the smoothing action of surface tension limits 
the size of perturbations that one can apply to the bubble 
shape. Surface tension must therefore be partially cir- 
cumvented in experiments that require larger azimuthal 
shape perturbations. 

Inflating the bubble in a short time, on the order of 
40 ms, prevents surface tension from smoothing out the 
longest-wavelength asymmetries in the bubble shape — 
those produced by the nozzle shape itself. Figure [5Jd 



shows the apparatus for producing these rapid "bursts" 
of gas: A reservoir is connected to the nozzle by a length 
of plastic tubing, which is constricted by a small solenoid 
pinch valve. The pinch valve was chosen to minimize the 
displacement of gas due to the valve's operation. The 
reservoir, whose volume may be adjusted from the ~1 mL 
volume of the tubing up to 16 mL, is charged with a 
pressure regulator to 0.3-5.1 kPa above the hydrostatic 
pressure at the nozzle (pressures were set relative to am- 
bient atmospheric pressure; values are quoted in this pa- 
per using a typical hydrostatic pressure at the nozzle of 
0.4 kPa). In the present work, the reservoir volume is 
16 mL unless otherwise specified. Meanwhile, a separate 
line supplies gas to the nozzle at a constant, low flow rate, 
inflating a bubble quasistatically. In the early stage of 
this slow inflation, a convex cap protrudes from the noz- 
zle, corresponding to increasing gas pressure as it grows. 
The burst is triggered when the output voltage of a pres- 
sure sensor (Honeywell DC002NDC4) reaches a pre-set 
threshold, corresponding to a particular interface shape. 
To create the burst, a delay generator momentarily opens 
the solenoid valve, triggers the high-speed cameras, and 
blocks re-triggering for ~20 s, long enough for the most 
violent liquid motion from the burst to subside. The ap- 
paratus thus satisfies both requirements for repeatablc 
bursting: a consistent pulse of gas delivered when the 
interface is at a consistent initial shape. In practice, con- 
secutive bursts are very repeatable. However, perhaps 
owing to details of pinch valve operation, quantitatively 
reproducing a pinch-off behavior between runs is difficult, 
especially for high-pressure bursts, indicating a sensitiv- 
ity to aspects of the experimental protocol or apparatus 
that have yet to be fully controlled. 



III. EFFECTS OF GAS DENSITY 

The gas inside the collapsing neck is not required for a 
model of the singularity at pinch-off [3, 0, EH 1 but prior 
work has found that it does affect the dynamics nonethe- 
less [H, EH EH EH • Furthermore, the role of gas in the 
production of the small satellite bubbles that remain af- 
ter pinch-off Q is unclear. Here I report on experiments 
that independently vary the composition of the gas inside 
the bubble and the ambient pressure of the experiment 
(thereby varying both gas density and static pressure in 
the liquid), in order to elucidate the role of gas as it 
changes the neck shape near the singularity, and leads to 
the formation of a satellite bubble. 

My experiments vary the composition and pressure of 
the gas, achieving gas densities from 0.039 to 5.5 g/L — 
thereby varying the density ratio p g /pi from 3 x 10~ 5 to 
4 x 10~ 3 , and the Atwood number (pi ~ p g )/(pi + p g ) 
from 1.000 to 0.989. (All gas densities are calculated 
from Ref. [22| for room temperature.) 

Measurements of the evolution of i? m i n for the different 
gases studied here confirm that changes to gas composi- 
tion and density do not strongly affect the inertial pinch- 
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FIG. 3: Minimum neck radius 7? m i n as a function of time 
until break-up r for axisymmetric pinch-offs of N2 (o, density 
1.14 g/L, 3 events), low-pressure He (x, 0.0411 g/L, 6 events), 
and SFe (+, 5.63 g/L, 7 events), t*, and thus r, is determined 
for each movie by fitting with Eq. [T] with a = 0.56 (dashed 
line). Only points below R m in ~ 120 (im were used for fitting. 
The consistency of these data show that the changing gas 
density does not strongly change the character of the singular 
dynamics over most of the collapse — however, SF6 points 
below 60 fj,m could not be fitted with an exponent of 0.56. 

off dynamics over most of the collapse, so that the effects 
of gas inertia may be considered as a perturbation to 
the symmetric void collapse dynamics. Figure [3] presents 
Rmin vs. t data for axisymmetric pinch-off of 3 gases from 
a 6 mm-diameter circular nozzle in water: dry N2, He, 
and SFq. The initial transient, where surface tension is 
important, persists down to i? m i n ~ 300 zxm; below that, 
pinch-off is an inertial collapse described by Eq. [TJ with 
a = 0.56 ± 0.03. This exponent was previously found to 
hold for quasistatic pinch-off from circular nozzles with 
diameters of 3-8 mm [6] . Notably, the SFg data in Fig. [3] 
differ from the data for other gases, in that points below 
Rmin ~ 60 fim are described poorly by a power law with 
a = 0.56, suggesting new dynamics at small i? m i n due to 
gas inertia Alternately, the effective power law of 

the entire collapse could be higher, in a range centered 
ona = 0.58. The anomaly at small r may also be cor- 
rected by adjusting the pinch-off time i», at the expense 
of introducing a significant deviation near r = 20 /its. 

A. Effect of gas density on neck shape and vertical 
motion 

Rather than being symmetric about its minimum, as 
one may expect from a leading-order analysis the 
collapsing neck of a gas bubble has a slight vertical 
asymmetry, with a minimum that moves upwards in 
time [ll|, Hil, [23J] . Burton and Taborek recently stud- 
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FIG. 4: A pressure drop in the upper portion of the bubble is 
responsible for the upward gas flow. The pressure in the gas 
supply line is plotted as a function of the time r before pinch- 
off is seen on video, for He at 0.17 g/L as it disconnects from 
a round, level, 6 mm-diameter nozzle. The pressure drops by 
9.0 Pa in the final 10 ms of pinch-off. Pressure is relative to 
the value measured 400 ms after pinch-off, once vibrations of 
the residual interface have died down. The data for SF6 are 
not significantly different, except for a lag due to the slower 
speed of sound. 



ied bubble pinch-off for gas-to-liquid density ratios from 
~10 -3 to 1, and found that this asymmetry and mo- 
tion became more pronounced for higher density ratios; 
they attributed these effects to a "jet" of gas that flows 
through the neck minimum 10]- Gekle et al. simulated 
and experimentally measured this same gas flow, and 
showed that during the collapse of a void created by the 
impact of a solid object on water, a slight gas pressure 
imbalance between the regions above and below the min- 
imum drives the upward gas flow [111 ]. Late in collapse, 
the gas reaches velocities much greater than those of the 
water and exerts a force on the liquid at the narrowest 
part of the neck, thus creating the upward motion and 
asymmetry. 

By placing a gas pressure sensor at the nozzle, I have 
determined a lower bound on the pressure drop that 
drives this gas flow through the neck minimum. Gas 
pressure measured immediately before pinch-off is shown 
in Fig. |U As expected, the pressure gradually decreases 
as the top portion of the bubble rises into water with 
lower hydrostatic pressure; this pressure drop acceler- 
ates significantly in the final 20 ms of pinch-off. As in 
the work of Gekle et al. [ll| , the present estimate of the 
pressure drop responsible for the gas flow is small com- 
pared with the absolute pressure in the neck — 1 part 
in 10 as measured here, and of the order of 1 part in 
10 2 in the simulations of the prior work [ll|. My mea- 
surement is a lower bound on the magnitude of the final 
pressure difference: because of the limited time resolu- 
tion of Fig. 01 and the location of the pressure sensor at 
the end of several cm of tubing, some of the rapid change 
in pressure just before the singularity might be masked. 
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Because this pressure drop is primarily due to the large, 
rising free surface above the neck minimum, and not to 
changes near the nozzle, it drives an upward gas flow 
through the neck minimum. 

Here I show that while this rapid gas flow may play a 
role in the final moments of pinch-off, it does not explain 
the vertical motion and asymmetry observed much earlier 
in the collapse. Instead, it is the gas's role in setting 
initial conditions, at the onset of the inertial collapse, 
that allows it to affect vertical asymmetry at these early 
times. 

To investigate the effect of gas density on the verti- 
cal asymmetry of the neck shape, I measure and com- 
pare the upper and lower cone half- angles of the neck 
profile, as indicated in Fig. [5£l The angles U and 9i 
are obtained by fitting a straight line to the regions of 
the profile for z between 7i? m i n and 12i? min away from 
the minimum. This measurement technique was chosen 
for its relative insensitivity to the details of the method, 
and the repeatability of its results; the measured angles 
would be smaller near the minimum, owing to R{z) be- 
ing nearly parabolic in that region. Fig. [SJd shows that 
under all conditions, the profile becomes more slender as 
Rmin decreases in the approach to pinch-off, consistent 
with earlier observations [14|, [l8[ . 

The asymmetry 6 U — 81 , which is apparent in the photo- 
graph of Fig. [5^, shows a strong dependence on gas den- 
sity, as shown in Fig. [5]:. Lowering gas density greatly 
reduces the difference between 6 U and primarily by 
making the upper shape more slender at pinch-off. Fig- 
ure [SU shows that gas density makes the greatest differ- 
ence when the neck is most constricted, at i? m i n < 20 /im 
and the velocity i? m i n > 7.5 m/s, as seen in the growth 
of asymmetry for the highest-density SF 6 (5.5 g/L) - 
whereas the less dense gases show a neutral or slightly 
negative trend of asymmetry in this late stage of collapse, 
consistent with recent simulations of pinch-off with ver- 
tical asymmetry and no inner fluid |24(. The qualitative 
difference for the highest-density gas at small i? m i n , as 
is also seen in the evolution of i? m i n (r) in Fig. [3J sug- 
gests that inertial effects change the dynamics late in the 
collapse. 

These results suggest that for most of the collapse, 
vertical asymmetry is due to a mechanism other than 
gas flow through the neck minimum, since there is an 
effect even at the much lower velocities earlier in the col- 
lapse. To explore this possibility, I turn to a related be- 
havior, the vertical motion of the neck minimum z m ; n , 
that is more easily measured at early times — includ- 
ing in the transient regime before inertial collapse, when 
surface tension is still important. Figure [6£i shows con- 
siderable, gas-dependent motion of z m i n before the onset 
of the inertial collapse at r ~ 1 ms (R m i n ~ 600 /im). 
Figure [Hb shows that this dependence carries over into 
late times. Around i? m i n ~ 50 /im (i? m i n ~ 3.6 m/s), the 
slope steepens further for pinch-off with the highest den- 
sity, consistent with the results for asymmetry in Fig.[5ji. 
These observations show that much of the vertical asym- 
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FIG. 5: (Color online) Influence of gas density on the ver- 
tical asymmetry of the neck, (a) Measurement of cone an- 
gles 6 U and For consistency and repeatability, cone angles 
are computed by fitting a straight line to the region of the 
profile between 7-R m i n and 127? m i n away from the minimum, 
(b) Lower cone angle 9i as a function of neck minimum ra- 
dius -Rmin, for 3 gas densities: lowest-vacuum SFe (5.5 g/L, 
circles), highest-vacuum SFg (1.5 g/L, squares), and highest- 
vacuum He (0.039 g/L, triangles). Time proceeds from right 
to left; thus the neck grows more slender during collapse, (c) 
Vertical asymmetry 6 U — 8i as a function of gas density p g , 
sampled at i? m i n = 10 /im (circles), 20 /im (diamonds), and 
60 /im (triangles). The data suggest a small residual asym- 
metry that is not due to gas inertia, (d) Vertical asymmetry 
as a function of neck minimum radius i? m i n for the 3 gas den- 
sities shown in (b). Higher gas density gives rise to a growth 
in asymmetry late in collapse. Earlier, asymmetry appears to 
evolve from an initial condition at large -R m i n that is deter- 
mined by the gas density, but with a trend that is independent 
of the gas. 



metry is present in the initial conditions of the inertial 
collapse. 

I also note that this vertical motion of z m i n causes an 
increase in the final height of the neck minimum. This 
effect leads to a problem at low gas densities: as the up- 
ward motion is nearly eliminated, a smaller volume of 
residual gas is left at the nozzle after pinch-off. Because 
this residual bubble vibrates after it is violently deformed 
by pinch-off, its flatter shape and lower height mean that 
the interface may depin and fall below the nozzle rim 
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the neck minimum. The vertical symmetry of the ini- 
tial conditions is presumably broken by a combination 
of buoyancy and geometry (the fixed nozzle below, and 
only a free surface above) which would be altered by a 
change in gas density. That this breaks the symmetry 
of the subsequent collapse may be understood by con- 
sidering the slenderness approximation discussed above, 
in which dynamics at different z are roughly indepen- 
dent. In this sense, the inner fluid acts as a perturbation 
that is remembered by the singular dynamics. By con- 
trast, the effect of gas flow through the neck observed for 
higher p g /pi and small i? m i n requires the introduction of 
new dynamics late in collapse, and so does not consti- 
tute a memory of initial conditions. This latter effect 
was prominent in the more energetic pinch-offs studied 
by Gekle et al. [llj . though that work did not investi- 
gate the role of the gas density. The work of Burton 
et al. 10], on a nozzle pinch-off experiment similar to 
the present one, varies the gas density in the transition 
from the perturbative regime to pronounced asymmetry 
and strongly altered dynamics. Those experiments only 
marginally overlap with the present work, which covers 
the weakly perturbative regime of extremely small p g /pi, 
and so can describe pinch-off and identify the role of gas 
in the near-absence of inertial effects. 



FIG. 6: (Color online) Upward motion of the neck minimum 
is driven by a vertical gas flow as a bubble disconnects from 
a round, level, 6 mm-diameter nozzle, (a) Height of the neck 
minimum above the nozzle as a function of minimum radius 
R mi „ for He at 0.029 g/L (•) and SF 6 at 5.7 g/L (o). In- 
creasing gas density causes an upward motion of the neck 
minimum, which for SF6 is more pronounced and begins at 
much larger Rmi n . r ranges from 12 to 0.1 ms over the course 
of this plot, (b) Height of the minimum as a function of i? m i n 
closer to pinch-off, showing a singular motion at late times 
that increases with gas density. z m i n is compared with its 
extrapolated value at Rmin = 0, z P mch, to make the curves' 
different slopes more apparent. From top to bottom, the 5 
curves correspond to He at 0.041 g/L, He at 0.16 g/L, SF@ 
at 1.6 g/L, N2 at 1.1 g/L, and SF6 at 5.6 g/L. The range of 
-Rmin plotted represents the final 100 ps of pinch-off. 



during the vibrations. Any depinning interferes with the 
repeatability of pinch-off (see Sec. Ill A| ). and is the practi- 
cal limitation on precision low-pressure experiments. For 
He gas, it causes difficulty at pressures well above the 
vapor pressure of water, even above 10 kPa. 

These results show that gas density has an observable 
effect on neck shape at all stages of collapse — even when 
the minimum neck radius is 2000 pm, and the maximum 
liquid velocity is just 0.05 m/s. I argue that this de- 
pendence is due to the effect of gas on initial conditions 
before inertial collapse, including the effect of transient 
dynamics dominated by surface tension. These initial 
conditions influence the subsequent evolution, and only 
at late times is the shape changed by gas flow through 



B. Effect of gas density on satellite bubbles 

These investigations also give insight into the possi- 
ble role that gas plays in the creation of satellite bub- 
bles. These bubbles, of diameter 10 pm or smaller, are 
observed near the point of disconnection after pinch- 
off [1, [H,[l^]. Their origin is not well- understood. The 
bubbles may contain vapor or gas that was released from 
the surrounding liquid by low dynamic pressure, as sug- 
gested by recent results of Burton et al. [l(| ■ Alternately, 
they may contain the same gas that comprises the main 
bubble, due to a mechanism that would bifurcate the 
neck minimum just before pinch-off, like that proposed by 
Gordillo and Fontelos @. Here gas, which flows rapidly 
in the vertical direction at the end of collapse, may pro- 
vide the necessary coupling in z at the small scales of 
satellite formation. In their model, when gas is forced 
out of the collapsing neck, there is a stagnation point at 
the minimum, and gas inertia lowers the dynamic pres- 
sure just above and below it, creating the bifurcation. 

Happily, these two possibilities — satellites from gas 
in the surrounding liquid, or from the motion of gas al- 
ready in the neck — produce two different dependencies 
on experimental conditions. In the former case, lowering 
the ambient pressure of the experiment would make the 
satellites larger, as more gas is drawn from the liquid; 
degassing the water would make them smaller. In the 
latter case, lowering the density of the gas would dimin- 
ish its ability to deform the neck shape, as discussed for 
the case of vertical symmetry above. Control over the 
ambient pressure, gas composition, and dissolved gases 
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makes these experiments possible. 

In the experiments using degassed water reported here, 
satellite bubbles were observed at all gas densities. There 
is a weak dependence of the satellite bubble size on gas 
density: the bubbles are ~10 (im diameter for SFg, and at 
most 3.5 /im diameter (the size of a pixel) for He. These 
experiments rule out dissolved gas as a source of satel- 
lite bubbles, both because the water was degassed, and 
because satellite bubbles are not larger at low pressures, 
when gas would come out of solution more readily. 

A previous effort by Burton et al. varying gas density 
instead observed no satellite bubbles at all, suggesting 
that the use of degassed water prevented satellite forma- 
tion [To| • This is not consistent with the results reported 
here. Although the mechanism developed by Gordillo 
and Fontelos [9[ shows how inertial effects of the gas may 
lead to the creation of a satellite bubble from gas already 
present in the neck, this particular mechanism posits a 
stagnation point at the neck minimum, inconsistent with 
the strong vertical gas flow discussed above. The results 
reported here suggest that satellite bubbles are due to 
the action of the gas already inside the neck, but there is 
presently no model consistent with these findings. 

C. Discussion 

The above results show that reducing the density of 
the inner fluid results in pinch-offs with smaller vertical 
asymmetry. However, over most of the collapse, the ver- 
tical asymmetry is due not to rapid gas flow through the 
neck — for which the case of zero gas density is easily un- 
derstood — but instead to initial conditions, where the 
effects of gas are weak and depend on the density differ- 
ence (pi — p g ). In particular, during quasistatic inflation, 
when buoyant and capillary forces are balanced, lengths 
in the system are non-dimensionalized by the capillary 
length 

lc = J , 7 , (2) 
V (Pi - P9)9 

where 7 is surface tension and g is the gravitational accel- 
eration. The analysis of Longuet-Higgins et al. [4J shows 
that the bubble shape at the end of inflation has a non- 
trivial dependence on the dimensionless nozzle size D/1 C1 
however, and so one may expect changes to the conditions 
at the end of inflation, the transient evolution, and the 
subsequent inertial collapse, beyond a simple rescaling. 

Fundamentally, the vertical symmetry of pinch-off is 
broken in two ways: gravity, which causes a hydrostatic 
pressure gradient, and geometry, specifically the presence 
of a fixed contact line at the bottom end of the bubble 
and a free surface at the top. I attempted to change 
the asymmetry by using an alternate geometry: an ap- 
paratus that replaces the free surface at the top of the 
bubble with a second, downward-pointing nozzle, so that 
the static bubble neck is created between the two paral- 
lel orifices. Pinch-off is initiated by slowly lowering the 



pressure inside the bubble. The use of this apparatus, 
however, alters the collapse behavior i? m i n (r), because 
surface tension holds the neck open until a much smaller 
^?min- This leads to poor observations of the singular dy- 
namics. In general, it is more difficult to isolate the con- 
tributions of gravity and geometry to symmetry-breaking 
than it is to identify the role of gas. 

IV. AZIMUTHAL MEMORY 

A key form of memory in bubble pinch-off is the preser- 
vation of azimuthal asymmetry present in the initial or 
boundary conditions. A perturbation to the initial con- 
ditions that breaks the azimuthal symmetry is remem- 
bered in the shape of the neck near pinch-off @ . A large 
class of these perturbations is addressed by the theory of 
Schmidt et al. [7], which considers azimuthal perturba- 
tions in a single horizontal slice of the neck, and describes 
the shape of that slice in terms of harmonic modes on a 
circle, 

r(0, t) = R(t) + J2 cos(<£„(i)) cos(n0 - 0°), (3) 
n 

where R(t) is the average radius of the cross-section, 
a n (t) is the amplitude of each mode, <f> n (t) is the phase 
of the oscillation, and 8 is the azimuthal coordinate. 0° 
selects the orientation of the perturbation, and is defined 
on [0, 2n/n). For a n -C R, these modes vibrate as 4> n (t) 
evolves during collapse, as illustrated for the n = 2 and 
71 = 3 harmonic modes in Fig. [7] Near the singularity, 
the phase evolves as 

4> n (t) = V^lHR{t)/R(t = 0)), (4) 

so that as R — > 0, in the limit of small a n , the shape 
goes through an infinite number of oscillations in finite 
time. Concurrently, the time derivative a n — > 0, so that 
any finite perturbation becomes more important as the 
collapse proceeds — the absolute amplitude of the per- 
turbation remains constant even as the average bubble 
radius shrinks to zero. 



A. Observations of the n — 2 mode 

Experiments performed by Schmidt et al. [7[ excited 
the n = 2 mode, and showed that the experimentally- 
observed vibrations agreed with the predictions of this 
theory. This same agreement is shown for a similar ex- 
periment in Fig. [8^. Here, bubbles are blown quasistati- 
cally from a slot-shaped nozzle, similar to that shown in 
Fig. [7]:. (For this particular experiment, I used a shorter 
4.6 mm-long version of that nozzle and found that in 
the quasistatic regime, results are indistinguishable from 
those of the longer slot.) From the two neck profiles mea- 
sured by the cameras in each frame of the movie (see ex- 
perimental methods, Sec. IIIBI) . the average radius R{z) 
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FIG. 7: (a) Horizontal neck cross-sections showing azimuthal 
perturbations [7J, and how they are measured. The length 
(A — B)/2, measured with two cameras, is the n — 2 asym- 
metry S2, and the average radius R = (A + B)/4. The length 
(C — D)/2 is S3, where C and D are measured from the ge- 
ometric center of the neck, marked by a closed dot; here, 
R = (C + D)/2. In practice, C and D are determined by 
assuming that the geometric center is stationary during the 
pinch-off, which is valid for the leveled nozzle used here [f|. 

(b) Cross-sections at a later time, showing the nature of the 
shape oscillations (not to scale). The midpoint of the n — 3 
shape, as seen by a camera looking toward the top of the page, 
is shifted to the right, even though its geometric center has 
not moved. The intermediate shape for each mode is a circle. 

(c) Nozzle orifice shapes used to generate these perturbations. 

(d) The 2-plate nozzle, used to excite an n = 2 oscillation by 
placing two transparent walls near a circular orifice. 



is computed by averaging the two half-widths. The n = 2 
asymmetry S2 is computed at the height z m - m where R 
is smallest (see Fig. [7]) . S2 is plotted here as a function 
of the minimum radius -R m in = R(z m i n ), which serves as 
a clock variable that tracks the time r = t* — t until the 
singularity. i? m i n has the advantage that it may be mea- 
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FIG. 8: Vibrational memory of n — 2 azimuthal perturba- 
tions. The independent variable R m - m is the average radius 
of the neck at its minimum; time proceeds from right to left. 
Theoretical curves (continuous lines) are computed accord- 
ing to Schmidt et al. 0, [25[ , taking the initial condition at 
ftnin = 128 and 156 /xm for (a) and (b), respectively. Repre- 
sentative error bars are shown on a fraction of the points, for 
clarity, (a) The n — 2 asymmetry 5*2 oscillates due to a slot- 
shaped nozzle (see Fig. [7|. S2 > corresponds to elongation 
perpendicular to that of the nozzle. 108 separate pinch-off 
events are plotted, (b) S2 oscillates due to parallel verti- 
cal plates on each side of the nozzle. 5*2 > corresponds to 
elongation parallel to the walls. 34 events are plotted. The 
frequency is the same in both experimental configurations, 
but the phase is slightly shifted and the transient behavior is 
significantly different. 



sured independently for each camera frame — permitting 
many events to be plotted together without modeling and 
fitting a for each one. Because R m m(T~) may be mod- 
eled as a power law (Eq. [IJ, a plot of these data on a 
logarithmic r axis would have a similar appearance. I 
have verified that experiments performed with a nozzle 
rotated by 45° show no signal in £2 as expected, since 
the n = 2 oscillation is only in orthogonal directions, as 
in Fig. El 

As expected, these azimuthal perturbations do not sig- 
nificantly alter the underlying axisymmetric collapse dy- 
namics. Figure [5] shows the evolution of R m m in the 
n = 2 oscillatory collapses of Fig. [5J along with the ax- 
isymmetric collapse of a nitrogen bubble from a 6 mm- 
diameter circular nozzle, as was also plotted in Fig. [3] 
The perturbed pinch-off is described well by Eq. Q] with 
a = 0.56 ± 0.03. (The slot-nozzle data have a smaller 
prefactor j3 than the axisymmetric case, likely due to 
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FIG. 9: Minimum neck radius R m i a as a function of time until 
break-up r for pinch-off from a 6 mm-diameter circular nozzle 
(o, 3 events); and minimum average radius R m i n of the n — 2 
oscillatory pinch-offs of Figs. [5^ (x, 12 events) and b (•, 4 
events). All plotted events used N2 gas. r was determined 
for each movie by fitting with the power law -R m i n oc r° 6 
(solid line). Only points below i? m m ~ 120 ^tm were used for 
fitting. 



demonstrate this principle. Bubbles were released from 
a 6 mm-diameter circular nozzle placed between large, 
parallel vertical plates spaced 10 mm apart, as pictured 
in Fig. [7J The plates are transparent, to permit imag- 
ing with 2 cameras at 90° angles as in the slot-shaped 
nozzle experiment. This configuration excites the n = 2 
vibrational mode with an amplitude of ~1 /im, with the 
same frequency but a different phase from that in Fig. [5^, 
- showing that the specific perturbation selects both 
the phase and amplitude of this linear oscillation. The 
plates' continual presence as a boundary condition that 
alters dynamics, rather than as an initial condition as for 
the shaped nozzle, is apparent in the absence of the steep 
early transient that is seen in the vibrations from shaped 
nozzles. Repeating the experiment with the plates re- 
moved restores azimuthal symmetry. 

The solid curves in Fig. [5] are theoretical n = 2 curves 
calculated with the method outlined in Ref. Q , with the 
additional refinement of a viscous term 25]. The only 
adjustable parameters are the empirical amplitude and 
phase of the oscillation. Each curve describes data only 
after the neck reaches some Rmin ~ 400 /im, after which 
the collapse is dominated by inertia, making the theory 
applicable. The agreement of the predictions with the 
data is excellent. 



the smaller effective nozzle size of the slot-shaped nozzle, 
whose smallest dimension is only 1.6 mm. The slot- nozzle 
data are closer to those from a 4 mm-diameter circular 
nozzle.) 

The data in Fig. [5k. do not oscillate around S2 = 0, but 
instead around ~1.5 /im. This offset does not appear in 
the data of Schmidt et al. 7} . The form of the data in the 
present work is the result of small refinements in analysis 
methods, which lead to excellent agreement among sev- 
eral runs of the slot-nozzle experiment, each comprising 
dozens of movies. The offset persists after exchanging 
the cameras and their optics, or when using a slot with a 
perpendicular orientation. This bias is likely due to the 
sensitivity of the vibrations to asymmetric external flows, 
and their corresponding pressure fields. Such flows are 
set up by the early evolution of the asymmetric bubble 
shape due to surface tension, before the inertial collapse. 

Figure [HJd directly demonstrates this sensitivity to 
asymmetric external pressure fields. I report for the first 
time that it is possible to excite azimuthal vibrational 
modes not only by altering the nozzle shape, but also by 
imposing asymmetric boundary conditions in the vicin- 
ity of pinch-off. Such boundary conditions cannot affect 
the pinch-off dynamics until the water is in motion, after 
which they create an anisotropic pressure field by im- 
peding water flow. This excitation is consistent with the 
theory of Schmidt et al., which holds that shape vibra- 
tions involve the pressure field of the surrounding water, 
extending out to some multiple of the initial neck ra- 
dius [7j. Figure [5b shows the results of an experiment to 



B. Observations of the n — 3 mode 

Here I report the first measurements of the second vi- 
brational mode, n = 3 [26]. A signal of this vibration 
is shown in Fig. 110a : this signal is measured according 
to Fig. [7] as an apparent side-to-side motion of the neck. 
Unlike the n — 2 data shown in Fig.[8j the n = 3 observa- 
tion comprises more than a full period of the oscillation. 
The principal development required for this observation 
was a means of generating an n = 3 perturbation large 
enough to be detected optically. The smoothing effect 
of surface tension, described above in relation to n = 2 
perturbations, is stronger at the shorter wavelength of 
n = 3. Experiments with bubbles inflated quasistatically 
from an n = 3-shaped nozzle (see Fig. [7J are thus un- 
able to excite that vibrational mode observably. Instead, 
I used small bursts of No g described above, to 

generate sufficiently large perturbations (in this case, a 
burst pressure of 7.5 kPa with no reservoir volume). De- 
spite the more rapid bubble inflation, the perturbation 
is very repeatable; Fig. [PPl comprises 21 separate events. 
No n = 2 signal is observed for an n = 3 vibration, and 
vice-versa, consistent with pure single-mode vibrations. 

Both the analysis of theory and the measurement tech- 
niques require a n -C R. This condition becomes rele- 
vant in the n = 3 oscillation of Fig. [TUb : at small R m i n , 
the relatively large oscillation amplitude (a 3 ~ 20 ^m) 
causes a disagreement between measurement and the- 
ory. Crucially, the successful measurement of both S3 
and _R m in requires a cross-section that is entirely convex, 
as in Fig.[JJ According to Eq.[31 this condition is satisfied 



11 




-1-lJ , I I '• I I I 

10 1 10 2 10 3 

Neck size (^m) 




200 /jm 



FIG. 10: Vibrational memory of an n — 3 azimuthal pertur- 
bation due to a 3-armed nozzle (see Fig. [7J). (a) Vibration 
signal 5*3 as a function of neck minimum radius -Rmin. S3 is 
subject to an arbitrary offset, as in Fig. [7] To better show 
the consistent phase and amplitude of the 21 events plotted, 
a different offset is fitted to each event. For -R m m > 50 fim 
(shaded vertical line), the neck is sufficiently large compared 
to the perturbation that the cross-section remains convex (see 
insets), and so 5*3 may be accurately measured. To the left 
of this line, the cross-section is concave and the measurement 
no longer corresponds to 53. This change in geometry is be- 
cause the perturbation amplitude remains nearly constant as 
the neck size decreases, (b) The valid portion of the S3 data, 
with a theoretical curve as in Fig. [S] The initial condition was 
taken at -Rmin = 280 /jm. (c) Frames from one of the movies 
plotted in this figure. Due to the n = 3 vibration, the ap- 
parent center of the neck oscillates from side to side in time. 
In the final frame, with -R m i n = 20 /im, the neck shape above 
and below the minimum recapitulates this prior side-to-side 
motion. 

so long as 

o 3 (t) cos(0 3 (t)) < (1 - V3/2)R min (t). (5) 

Using the parameters of the oscillation in Fig. [TUb . this 
predicts a marked transition to concavity at R m in ~ 
50 fim, as illustrated by the insets of Fig. [TTJk . This cut- 



off is used to replot the data, with the oscillation curve 
predicted by theory, in Fig. [TUb . The theoretical anal- 
ysis, performed in the same manner as for the n = 2 
oscillation, is in excellent agreement with the valid range 
of experimental data. The cutoff also means that a mea- 
surement of R ni i n (t) has a small dynamic range and con- 
strains the power law poorly, but it is consistent with 
Eq.CD for a = 0.56. 

Remarkably, the n = 3 vibration of Fig. [10] is also 
apparent in still images of the collapsing neck. This is 
because of the nearly 2-dimensional dynamics responsible 
for the shape oscillations, which permit the vibrations to 
be modeled in a single horizontal slice, independently of 
those at other heights 0] ■ Because the initial neck shape 
perturbation varies slowly with vertical position z, slices 
of the neck at different heights are executing nearly the 
same vibrations but at different times — determined by 
their local average radius R. I have found that for all of 
the regimes of vibrating pinch-off presented in this paper, 
the oscillation phase and amplitude vary by only ~10% in 
the region \z— z m j n | < 200 /im. Thus slices above or below 
the neck minimum, having a size R > -R m ; n , recapitulate 
the oscillations already performed at the minimum. The 
last frame of Fig. [TUb shows such a neck shape, in which 
the apparent n = 3 side-to-side wobbling in time is mani- 
fest as a "kinked" neck. This kind of neck deformation in 
discrete steps may account for some of the neck shapes 
observed with other burst perturbations, as illustrated 
in Fig. 4a of Ref. @. My experiments for n = 3 have 
isolated this striking behavior, which should be a generic 
feature of vibrational memory with odd n. 



C. Nozzle tilt perturbations 

Notably, azimuthal shape vibrations are not excited 
by the tilting perturbations of quasistatic pinch-off, first 
described by Keim et al. Such a perturbation is cre- 
ated by tilting the nozzle orifice by an arbitrarily small 
angle from the horizontal; an effect for a 6 mm circu- 
lar nozzle can be observed for tilts as small as 0.05°. 
Figure [TTk shows that tilting the nozzle causes the neck 
minimum to translate horizontally away from the direc- 
tion of tilt during collapse. The final neck shape has a 
correspondingly lopsided appearance. At pinch-off, the 
neck minimum is thinner in the asymmetric view than in 
the symmetric one, reminiscent of the cross-sections pro- 
duced by an n = 2 oscillation. But while this perturba- 
tion is an excellent example of the memory of azimuthal 
asymmetry in pinch-off, I do not observe any associated 
oscillatory motion. The tilting memory is also different 
from the oscillations in that it involves the entire three- 
dimensional flow around the neck: the neck's horizontal 
motion has a strong additional dependence on z, meaning 
that the perturbed dynamics cannot be approximated as 
two-dimensional. By contrast, the shape oscillations near 
the minimum depend only on the local R(t, z), resulting 
in the behavior described above for n = 3, wherein the 
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FIG. 11: (a) Horizontal translation of the neck minimum 
as a function of neck Siz6 Rmin for a 6 mm-diameter noz- 
zle leveled to within 0.02° (+), and tilted by 0.1° (A) and 
0.75° (o) (4 movies each). Motion is away from the direc- 
tion of tilt. Unlike that oscillation, translation from tilting is 
monotonic. (b) By contrast, the signal from placing a single 
vertical wall 8 mm from the nozzle center (•, 14 movies) is 
weakly non-monotonic, changing at the onset of inertial col- 
lapse (~400 /jm). Here, motion is initially toward the wall. 



neck above and below the minimum repeats the oscilla- 
tions recently performed at the minimum. 

The tilting perturbation may be compared with the 
case of placing a single vertical wall ^8 mm from the 
nozzle center — a variant of the 2-plate geometry shown 
in Fig. [7J1. As shown in Fig. ITTb . this perturbation leads 
to a pinch-off with a weakly non-monotonic side-to-side 
motion. Unlike with the tilting perturbation, the motion 
depends primarily on R(r,z). In this way, it is like the 
n = 2 and n = 3 oscillations described above. At this 
time, neither a model of the tilting perturbation nor of 
the single-wall perturbation has been developed. 
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FIG. 12: (Color online) Varying burst pressure varies the os- 
cillation parameters. Oscillation signal S2 is plotted vs. aver- 
age neck radius at the minimum, Rmi n . Each curve represents 
a burst pressure, labeled in units of kPa. Here Rmin has a lin- 
ear scale, because of the importance of behavior before the 
inertially-dominated power-law collapse. At lower pressures, 
each increase in the pressure advances the oscillation phase 
and increases its amplitude; this effect becomes more com- 
plicated, and oscillations less sinusoidal, at higher pressures. 
As in Fig. [HJ S2 > corresponds to elongation perpendic- 
ular to that of the nozzle. In the lower left, the oscillation 
for quasistatic inflation (i.e. zero burst pressure) from Fig.[HJi 
is plotted with its S2 multiplied by —10 and offset for com- 
parison. Increasing the burst pressure gradually changes the 
phase of the oscillation, and generally increases its amplitude. 
For each pressure, data points are plotted with a smoothed 
curve to aid the eye. The points for 0.35 kPa comprise 15 
pinch-off events; each other curve comprises 2. 



D. Vibrations and outcomes of large perturbations 

Schmidt et al. find that the amplitude of an oscillatory 
perturbation, as an absolute length scale, becomes fixed 
as the overall size of the neck collapses to zero 0] ■ Thus 
any finite perturbation will eventually dominate the neck 
shape. This regime has been analyzed by Turitsyn et 
al. [8J , who show that the cross-sections in this regime are 
not described by Eq. [3l which applies only when S n -C R 
as was the case in Fig. [JJ 

The bursting apparatus can create perturbations large 
enough to make this behavior visible, and allows access to 
a variety of oscillation parameters and outcomes by ad- 
justing the pressure of the burst. Figure IT21 demonstrates 
this control by showing the n = 2 oscillations from a se- 
ries of movies taken with a range of burst pressures. At 



low pressures, increasing the pressure results in larger 
amplitudes and advances the oscillation phase. While 
the data are not recognizable as memory oscillations at 
higher pressures, the low-pressure bursts still oscillate as 
though they are in the limit of small perturbations: the 
data for a burst of 0.35 kPa are remarkably similar to 
the quasistatic data that are also plotted in Fig. even 
though the perturbation amplitude is larger by an order 
of magnitude. 

It is only when this 0.35 kPa pinch-off has collapsed 
down to the scale of the perturbation that there is a 
dramatic difference, as shown in Fig. (T3J Instead of the 
pinch-off appearing to end at a single point, two sides of 
the neck make contact and merge, as in the coalescence of 
two liquid drops [|| . The perturbation scale thus cuts off 
the singularity, replacing simple disconnection at a point 
with a different topological change [fUl]. The particular 




FIG. 13: Outcome of a small burst of N 2 (0.35 kPa) from 
a slot-shaped nozzle, as seen simultaneously from a direction 
perpendicular to the slot (left) and parallel to it (right). In 
the left image, there is a satellite bubble in the center of the 
hole. Scale bar is 0.5 mm. 



structure and dynamics of the coalescence outcome are 
detailed in Sec. IIV El below. Large perturbations alter 
the singularity in every case, though I do not always ob- 
serve a coalescence — for some conditions, coalescence 
may be too small or fast to be imaged, or pinch-off may 
end with entirely different geometries, such as those sug- 
gested by recent simulations |27j . The present work, how- 
ever, focuses on the case of coalescence, because it cre- 
ates features in camera images {e.g. the appearance of 
a hole) that can easily be related to the topology and 
three-dimensional shape of the neck. 

By varying the oscillation parameters, changing the 
burst pressure also alters the way in which the singular- 
ity is disrupted. As the neck executes an n — 2 oscil- 
lation, the long axis of a cross-section changes orienta- 
tion by 90° for each half-period of the oscillation. Since 
increasing the perturbation amplitude will cause it to 
disrupt the dynamics sooner in the collapse, one may 
expect that amplitude to select the orientation of coales- 
cence events — an idea more fully developed by Turitsyn 
et al. @. Indeed, in these experiments, further increas- 
ing the strength of the burst results in coalescences at 
larger scales, but also changes the orientation of the coa- 
lescence event: while the moderate-sized event of Fig. [13] 
takes place along a plane parallel to the slot nozzle, the 
larger event ofRPfk has a similar geometry but in the per- 
pendicular direction. These observations indicate that 
memory oscillations happen even for large perturbations 
and -R m in that are not strictly within the inertial collapse 
regime of e.g. Fig. |S]that Schmidt et al. analyzed Q. 



E. Coalescence dynamics 

Examining larger (and therefore slower) coalescence 
events gives the clearest view of this outcome of strongly 
perturbed pinch-off. Figure [TJ] indicates that smooth, 
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flat interfaces from opposite sides of the neck are brought 
into proximity, momentarily forming a thin sheet of gas. 
Newton's Rings, seen in Fig. [Tib , and the estimated vol- 
ume of the resulting satellite bubbles, indicate that this 
sheet is on the order of 1 /im thick. This smooth contact 
is predicted by the analysis of Turitsyn et al. [8j] for large 
n — 2 distortions, as in the cross-section in Fig. [bib . 

Coalescence corresponds to the breakup of the thin 
sheet of gas. Instead of observing the process to proceed 
radially outward from the putative axis of symmetry, as 
it would for two liquid drops brought into contact |28l - 
l3lj . one can see one of two related behaviors: For the 
larger perturbations, as in Fig. Q3] a front of disinte- 
gration sweeps across the thin sheet, leaving hundreds 
of satellite bubbles. This behavior is reminiscent of the 
analysis of inviscid two-drop coalescence by Duchemin 
et al., who find a rapid sequence of breakups as the air 
sheet retracts [30j . For smaller perturbations, however, 
as in Fig. 1151 the sheet simply retracts around its perime- 
ter, leaving a central piece of the sheet that typically 
condenses into a single satellite bubble. Coalescence pro- 
ceeds rapidly, and with the exception of large events as in 
Figs. 1141 is not observed in more than one frame of video, 
suggesting a timescale of approximately 5 fis or less for 
events such as Figs. [T3] and [15] The implied rewetting 
velocity at which the front of breakups moves across the 
sheet is several m/s, consistent with balancing liquid in- 
ertia with the action of surface tension at the edge of the 
1 /im-thick sheet [3(j. Finally, an instance of another 
remarkable feature of large perturbations is the waviness 
of the outside neck edges in Fig. [15] This shows the neck 
shape's recapitulation of the large- amplitude n = 2 oscil- 
lations at the minimum, analogous to Fig. llOb . These co- 
alescence outcomes illustrate the rich three-dimensional 
structures that may result from simple perturbations at 
the beginning of pinch-off. 



V. CONCLUSIONS 

In this paper, I have presented experiments that fur- 
ther reveal some of the possibilities for perturbing the 
fully non-universal fluid pinch-off of gas bubbles under- 
water. Shape vibrations may be excited not only by 
changing the nozzle shape, but by breaking the azimuthal 
symmetry of nearby boundary conditions. The model of 
Schmidt et al. [7\ applies beyond the n = 2 mode previ- 
ously studied, agreeing with my observations of the n = 3 
vibrations that it predicted. When the perturbation scale 
is increased to a size that can be imaged more easily, these 
same dynamics are seen to produce fully 3-dimensional 
shapes, which require an alternative imaging technique 
such as X-rays [8|, [27| . These outcomes — including the 
remarkable coalescence dynamics that can leave hundreds 
of satellite bubbles — show not only the ability of this 
type of implosion singularity to retain information, but 
its ability to create complex shapes over a range of length 
scales from a simple perturbation on a single length scale. 
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FIG. 14: (a) The morphology of Fig. [13] is explained by observing a much stronger burst (5.5 kPa) with the same orthogonal 
cameras, producing simultaneous images "1" and "2." A similar geometry is formed but on a much larger scale, but note 
that the geometry has changed orientation. Scale bar is 0.5 mm. (b) Enlarged and contrast-enhanced detail from (a) shows 
that the "hole" is actually a ~1 /^m-thick sheet of gas, as indicated by the presence of Newton's Rings. Right: 87 jj,s later, 
breakup of the thin sheet proceeds from left to right, leaving hundreds of satellite bubbles, (c) Horizontal cross-section of the 
experiment, showing placement of the cameras used in part (a), and a neck cross-section predicted by Turitsyn et al. @] that 
exhibits smooth contact between opposite sides of the neck, consistent with formation of the thin sheet. 



My experiments further implicate inertial gas flow in the 
creation of satellite bubbles, but since the gas flow rules 
out the stagnation point that is required by the theory 
of Gordillo and Fontelos [9], the underlying mechanism 
by which gas bifurcates the neck minimum remains an 
open question. An understanding of satellite bubble for- 
mation might be of particular importance to applications 
in which one size of bubble is needed exclusively. 

The results presented here are representative of the 
richness of bubble pinch-off that has been revealed 
through the recent work of many investigators. This com- 
plexity underscores the potential for further insights into 
this remarkable class of singularities that retain informa- 
tion. 
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FIG. 15: 300 ns arc-flash exposures taken at different times 
in the course of two similar bursts from a slot nozzle. Scale 
bar is 0.5 mm. Left: Newton's Rings show a thin sheet as in 
Fig. 1141 Right: coalescence proceeds around the perimeter of 
the sheet, showing a process for producing a single, central 
satellite bubble as in Fig. 1131 rather than many bubbles, as 
in Fig. [H 

The gas within the bubble is likewise a remembered 
perturbation, in that it contributes, by way of initial con- 
ditions and transient dynamics, to the vertical asymme- 
try and position of pinch-off. With sufficient gas density, 
the effects of rapid gas flow through the neck minimum 
are also observed as new dynamics late in the collapse. 
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